Norm breaking in a queue — athermal phase transition 
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Standard order-disorder phase transition in the Ising model is described in terms of rates of 
processes of spin flips. This formulation aUows to extend numerous results on phase transition for 
sciences other than physics of magnetism. We apply it to the problem of norm breaking. Two 
strategies: to cooperate or to defect, form an analogy to two spin orientations. An example is a 
crowd at an exit, where to defect means to push others, and to cooperate is to move slowly together. 
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I. INTRODUCTION 

Phase transitions are ubiquitous in nature. The ferro- 
paramagnetic phase transition in Ising system- is an ex- 
ample cited perhaps most frequently. The concept of 
phase transition is valid and useful also in social sci- 
ences, but the example taken from magnetism seems 
to dominate minds. As an effect, we read on social 
temperature^i"— and social energy^, although these phys- 
ical quantities have no direct counterparts in sociology. 
Here we argue that in social and presumably also other 
sciences, it is more natural to describe a phase transition 
in terms of rates of processes, which appear in the prob- 
lem under consideration. Further, the standard magnetic 
phase transition can be also described in terms of rates 
of processes: in the simplest one-spin mean field model, 
there are two processes: the spin fiip from up to down 
and from down to up. This indicates, that the approach 
advocated here is more general. A similar scheme has 
been proposed by Galam and coworkers on iterative equa- 
tions, with discrete time^i^. Here we work with Monte 
Carlo simulations, supplemented with differential Master 
Equations. 

We do not claim that concepts of thermodynamics is 
not applicable in social sciences. In particular, consider 
a large ensemble of similar and independent events, each 
described with a variable x. Analysis of the probabil- 
ity distribution of x, when combined with the maximum 
entropy principle, can give a valuable insight into the na- 
ture of underlying processes. A recent example of this 
procedure can be found in Ref. To accomplish this, 
however, we need an assumption on the equilibrium state 
of the system, what is rarely firm in sociology and can 
only be justified ex post, from pertinence of obtained re- 
sults. On the other hand, rates of changes of investigated 
variables allow to reproduce time dependences of inves- 
tigated variables, and not only their equilibrium state. 

As an example, we chose the problem of norm break- 
ing at a crowded exit. Crowd simulation become an inde- 
pendent branch of computational science, with numerous 
applicationsiSrJ^. For our purposes, the system can be 
reduced to a set of cells occupied by pedestrians. We 



assume that the crowd is temporarily jammed, therefore 
pedestrians cannot move. However, each pedestrian can 
push other pedestrian before him, to activate the mo- 
tion as soon as possible and to reach the exit, or at least 
to force people before her/him to push as well. This 
small psychological advantage, however, can have fatal 
consequences, as the force exerted by particular pedes- 
trians accumulates at the exit and people there can be 
badly injured, if the crowd is large and the selfish pushing 
strategy prevails. This aspects allows to interpret the sit- 
uation as an example of the Prisonner's Dilemma—. The 
connection between crowd dynamics and game theory has 
been emphasized by Dirk Helbingi^. In particular, we 
expect that the phase transition in minds of pedestri- 
ans, which we investigate here, is coupled to the phase 
transitions discussed in theory of crowd dynamics, where 
the stop-and-go phase and crowd turbulence have been 
distinguishe d "'^^1-'^^ . 

The text is constructed as follows. In the next section, 
the Ising phase transition is described in terms of rates 
of spin-fiips. We provide an example, where the rates are 
taken from the standard thermodynamic, temperature- 
dependent formulation. For this case computer simula- 
tion reproduces the exact Onsager solutioi*ii. We pro- 
pose also a mean field formulation of the problem in 
terms of Master Equations. As a rule, the mean field 
approach gives a value of the Curie temperature, which 
is larger than the approximately true value obtained in 
the simulation. In the thermodynamic case, the mean- 
field solution is obtained which is slightly better than 
the Bragg- Williams solution but worse than the Bethe- 
Peierls solution^. Here, the criterion is the distance be- 
tween the obtained critical point and its exact value. Sec- 
tion mil is devoted to the crowd problem, as remarked 
above. As the full formulation of the problem on the 
square lattice contains many parameters, analytical cal- 
culations are presented for two one-dimensional cases, 
where i) two strategies are symmetric, ii) the homoge- 
neous states are absorbing. Whilst the fork bifurcation- 
is found in the case i), in the case ii) the bifurcation is 
transcritical. Numerical calculations are presented also 
for pedestrians in cells of a square lattice, with model 
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set of rates. In the same section we demonstrate that 
the relaxation time r goes to infinity at the critical sub- 
space of the spin-flip rates. Last section ITVl is devoted to 
conclusions. 



II. ATHERMAL ISING PHASE TRANSITION 

In physics, the ferro-paramagnetic transition of Ising 
spins s = ±1 at the square lattice is referred to most 
often. In this transition, the energy of the exchange in- 
teraction between nearest neighbours competes with the 
thermal noise. As a result, the ferromagnetic phase ap- 
pears at low temperature, whilst at high temperature the 
directions of spins are disordered. The same content can 
be expressed as follows. Spins flip from up to down and 
the oppositely, but the rates of these flippings depend 
on the states of spins in their direct neighborhood. In 
square lattice a spin up surrounded by four spins up flips 
rather rarely. On the other hand, a spin up surrounded 
by four spins down flips likely. Let us denote these rates 
as Wn, where n = 0, - ■ ■ ,4 marks the number of spins at 
nearest four nodes in the same state as the spin under 
consideration. Then, a spin flips from up to down with 
rate w„ if it is surrounded by n spins up and (4 — n) 
spins down. From the up-down symmetry, a spin down 
surrounded by n spins down flips to the up orientation 
with the same rate w„. 

In statistical mechanics, spins flip with the rates which 



fulfill the condition of detailed balance. For a spin up in 
the square lattice with n neighbours up, this condition is 
aWn — 6i/74_„, where a = (l + m)/2 is the probability of 
spin up, b = (1— m)/2 is the probability of spin down, and 
TO — Si/N , where summation is over N sites i. As the 
ratio a/b — '^^(""2) equilibrium (where 7 = exp(/3J), 
/? = 1 /T is the inverse temperature and J is the exchange 
integral), and the rates Wn do not depend on the fact if 
the system is in equilibrium or not, these rates should also 
fulfill the condition u'4_„ /wn = 74("-2). This is fulfilled 
in particular for the Metropolis algorithmic, where wq — 
wi — 1,W2 — 1/2, W3 = exp(— 4/3J) and W4 — When 
these rates are introduced to the Monte Carlo code, the 
obtained magnetization vanishes at W3 = 0.17, as shown 
in Fig. [TJ This gives the critical temperature Tp/J ~ 
2.27 , what agrees well with Onsager solution^^. 

The conditions of detailed balance fix the relations 
between w„ and Wz-n only, where z is the number of 
nearest neighbours. Then the relation between, say, 
and wi remains unspecified. In particular, we can set 
Wn — 7"*^^" However, in numerical calculations we are 
interested in high rates, to speed up the calculations, 
and therefore the Metropolis version is more convenient. 
In our text, the Metropolis rates are used for the sake 
of comparison of our results with the special case, where 
the detailed balance conditions are fulfilled. In general, 
the thermodynamic balance conditions are not valid, and 
our rates do not fulfill them. 
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The related mean field theory can be formulated in terms of the Master Equations^ 

§ = E „ + bwz^^). (1) 

at ^-^ n\(z — n]\ 

n=0 ^ ' 

To check its result again in the thermodynamic case, we can put Wn — 7^"'^". Substituting this to Eq. ([1]), after 
simple manipulations we get 

^ = _a(a/7 + 67)-+6(a7 + 6/7)^ (2) 

In this symmetric case, the stationary solution a ~b (m — 0) always exists, but it ceases stability when 

d da{m) 



dm dt 



= (3) 



at TO = 0, what gives tanh(/3J) = 1/z at Tc- This result is slightly better than the simplest Bragg- Williams mean field, 
where no correlations exist and zfiJ = 1. However, it is worse than the Bethe-Peierls result that tanh(/3J) = l/(z— 1) 
at Tc- This is because the states of next neighbours are not correlated with each other; in Eq. ([1]), their weights are 
given by the Bernoulli distribution. For z = 4, the stationary solution for the ordered phase is 

7-5 + 4t/ + 6y2 - 4y3 - + 2^{y - l)^{y + 1)2(5 + 2y + ^ 
' 

where y — ^"^ — exp(2/3J). 
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We could also try to use the 'Metropolis rates' wq = wi = 1, W2 = 1/2, W3 = x and W4 = a;^, to construct 
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m(0) exp(At). As r is usually defined by the relation 
m{t) — to(0) exp(— t/r), we have t = — 1/A. 

To conclude this section, Eq. JT]) provides the gener- 
alization of the mean field theory for non-physical (non- 
thermodynamic) cases. This means that the rates Wn do 
not fulfil the thermodynamic conditions of detailed bal- 
ance, but they are mutually independent. Also, their use 
does not imply any form of energy or temperature. As we 
have shown, phase transitions appear as a consequence of 
interplay of flows of probability between different states. 



III. NORM BREAKING IN CROWD 



FIG. 1: Mean value m of the variable s — ±1 (points) against 
parameter a; in a square lattice. When the rates are selected 
as Wo = wi = 1, W2 = 1/2, ws = X and W4 — x^, the 
system dynamics is equivalent to the Metropolis algorithm, 
where x — exp(— We obtain the curve of crosses (x) 
and the transition at s = 0.17. For the rates wq = 1, wi — 
2x, W2 — X, = x/2 and W4 = x/A we get the curve of 
pluses (+), where the ordered phase ends at about x = 0.06. 
The simulations were carried out for square lattice with A'^ — 
100^ spins and heUcal boundary conditions. The spins were 
visited in type-writer order. For given x parameter Tmax = 
10*^ Monte Carlo steps per spin were performed. The mean 
value of magnetization m{x) was evaluated as time average 
over the last Tmax/S = 5 • 10® steps. 



the Master Equation. In this case, the critical value of 
X is about 0.572. Substituting x — exp(— 4/3J), we get 
Tc = 7. 14 J. This is much worse even than the Bragg- 
Williams solution Tc = 4, while the rates w„ = 7^"^" 
give Tc = 3.91, the Bethe-Peierls solution is Tc = 2.89 
and the exact solution is Tc = 2.27 in units of J. This 
result is not strange, because in fact the detailed bal- 
ance conditions are not fulfilled in the stationary state 
obtained within the mean field theory. Instead, the con- 
tributing processes on the right side on Eq. ([T|) mutually 
cancel at the stationary state. When we put wi = 1, 
W2 = 1/2 and W3 = x, the processes of flipping of spins 
with one, two or three neighbours in the same state are 
much faster, than they should be if wq =^ 1, and the 
balance of processes is disturbed. 

The time t of relaxation to a stable state can be eval- 
uated by linearization procedure. As long as the state 
TO = is stable, the equation of motion can be written 
for small to as 



da{m) 



dt 



Xm, 



(5) 



where 
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d da{m) 



dm dt 

at to = 0, and A < 0. The approximate solution is m(t) 



(6) 



There is a countless number of theories of human 
behaviour, each one bearing some element of rcduc- 
tionism. Our choice is the Scheff theory of pride and 
shame^^'^^. There are two loops of behaviour, the shame 
loop and the pride loop, which relate to bad and good 
self-evaluation^. Accordingly, agents in the pride loop 
are more willing to select altruistic strategies, while those 
in the shame loop are likely to select the selfish one. 
Additionally, in the pride loop the communication with 
other agents is more efficient. Being in a loop, an agent 
prefers to remain there, hence the term; still, there is 
some probability to switch. In particular, a friendly com- 
munication with other agents can trigger the switch from 
the shame loop to the pride loop. 

When we project this theory to the problem of norm 
breaking, and to our example of the crowd behaviour 
in particular, it is clear that both the selfish strategy 
and the cooperative strategy have the property of self- 
stabilization. Further, the selfish strategy can spread; 
being pushed by other agent reduces the pride and drives 
the self-evaluation down. On the other hand, an example 
of a calm behaviour can bear an attitude to copy this 
strategy as well. This mechanism works better, if more 
neighbours behave in the same way---*2i. Concluding, it is 
justified to introduce the rates of strategy changes, which 
depend on the number of neighbours which play a given 
strategy. 

In principle, the process can be investigated in any 
space dimension, but the number of parameters Wn in- 
creases with the number z of nearest neighbours. To 
keep the analytical description simple, let us consider 
a one-dimensional queue, where a pedestrian has only 
two neighbours. As it is known, the ferro-paramagnetic 
phase transition in the one-dimensional Ising chain ap- 
pears only at T=0 K. Here its presence is a consequence 
of the mean field approach. Now two orientations of spin 
mean two options: the selfish strategy with probability a 
and the cooperative strategy with probability h — \ ^ a. 
In a linear queue, there are two neighbours only. In the 
'additive' case, i.e. if neighbours before and behind a 
given pedestrian are equivalent, the equation can be writ- 
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ten as 

— — c?[—aw2+hwa)^2ah(^—aw\+hw\)^\?[—awQ^-hw2)^ 

(7) 

where three consecutive terms on r.h.s. refer to three 
configurations of the neighborhood: two selfish neigh- 
bours (a^), one selfish and one quiet (2a&) and two quiet 
ones (6^). As we have described the change of each strat- 
egy with the same parameters, again the solution a — b 
is always possible. As before, the critical point can be 
identified from the instability of this solution. Denoting 
the ratio {wq — 2wi)/w2 as cj), we get the solution for m 
different from zero and stable if (A > 3. There 



(8) 
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otherwise m = 0. The counterpart of the critical temper- 
ature is then the value = 3. This is an example of the 
fork bifurcation. To have an ordered state, where most 
pedestrians apply the selfish strategy or most pedestrians 
apply the quiet (cooperative) strategy, we have to keep 
the rate wg large enough. This means that a pedestrian 
should change its state quickly enough when it is different 
than the state of both her/his neighbours. 
In our case, we get the relaxation time 



(9) 



W2(3 - (/)) 



and T tends to infinity as (p tends to its critical value </) = 
3 from below. This means, that the mean-field version of 
'critical slowing down' is present here. 

In the case of queue, we would like to consider also 
the case when a pedestrian makes a difference between 
the states ab and ba; it is perhaps not the same situation 
to be pushed from behind or just to observe that the 
neighbour in front of you is pushing somebody else. Also, 
the assumption on symmetry should be released; to start 
pushing and to stop pushing should be described with 
different rates. If these differences are introduced, the 
equation of motion is 

— = a^{-awaa + bUaa) + ab{-aWab + bUab) + 
ba{-awba + buba) + b^{-awbb + bubb), 

where Uxy {wxy) is the rate of start pushing (stop push- 
ing) when a pedestrian behind applies the strategy x and 
the pedestrian before applies y. When the equation PH)) 
is rewritten as 

^ = -a^Waa + b^Ubb - a^b(Wab + Wba ^ "aa) + 

ab^{Uab + Uba - Wbb) 

the number of parameters falls to four. In a generic case, 
neither m = 0, nor m = ±1 is a solution. The sta- 
tionary state can be found from the third-order algebraic 



equation by the Cardano method. As the number of real 
solutions can be either one or three, we expect that the 
related bifurcation is again of fork type. 

As a special case, let us consider the case when there 
is no outflow from the states aa and bb. Then, Waa = 
Ubb = 0, and the equation of motion reduces to 

^=2f^^'-(l-m'){D-C~miD + C)), (10) 

where C = Wab + Wba - Uaa and D ^ Uab + Uba - Wbb- 

There are three fixed points here, m± = ±1 and toq = 
{D — C)/{D + C). li D + C is positive, toq exists and 
is stable if both C and D are positive; having changed 
the sign of both C and Z?, we have the same rule. In two 
remaining areas either m_|_ or m_ is stable; the respective 
conditions are C < and D <Q. There is no more than 
one stable solution for a given set of parameters. The 
bifurcations are transcritical and not of fork type. As in 
the previous case, the linearization procedure allows to 
obtain the relaxation times analytically; for m+ (stable if 
C < 0) we get T = -1/C, for m_ (stable if D < 0) we get 
T = -l/D, and for mo (stable iiDC/{D + C) > 0) we get 
T — {D + C) / (DC) . In all cases, the relaxation time goes 
to infinity when the fixed point goes to the boundary of 
the area of its stability. Here these boundaries are C = 
and D = 0. We expect this behaviour of r also in more 
general cases. 

The order-disorder phase transition in one-dimensional 
systems is usually an artifact of mean field theory. On 
the other hand, a natural medium to simulate processes 
in crowd is a two-dimensional lattice. As we have seen, 
for the case of two symmetric strategies we need five pa- 
rameters Wn, n = 0, ■ ■ ■ ,4. In the presence of an ordered 
phase the rate w„ should decrease with n. This means, 
that there is a positive correlation between actions of 
neighboring pedestrians. Having this in mind, we have 
calculated numerically the magnetization on a square lat- 
tice, with the rates as follows: wq = 1, wi = 2x, W2 = x, 
W3 — x/2 and W4 — x/A. The choice of = 1 is equiva- 
lent to an assumption, that the strategy opposite to what 
all neighbours do is abandoned immediately. Other rates 
are set as that an increase of the number of neighbors who 
perform the same strategy by one leads to doubling the 
mean time when the strategy is continued. As we have 
no experimental data on these rates, the choice remains 
arbitrary. However, the decrease of Wn with n has some 
support in the results of observations of Mark Granovet- 
ter on the so-called critical mass effec t^^i^^ . The results 
of the simulation is shown in Fig. 1 as the curve of pluses. 
The approximate value of the transition point xc is 0.06. 
Note however, that for the latter red curve the parameter 
X has no thermodynamic interpretation. The mean field 
equation ((T|) for these rates gives the stationary solution 



m — a 



'4 + 5a; - 2^16 - 56a; -t- 85a;2 
15a; -4 ' 



(11) 



with the critical point close to a; = 0.571. Here again, its 
value is overestimated by the mean field model. 
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IV. CONCLUSIONS 

We have proposed a scheme of translation of thermal 
phase transitions to phase transitions related to an in- 
terplay of rates of different processes. This enables wide 
and easy application of results obtained in frames of sta- 
tistical mechanics to other sciences, where temperature is 
not a handy concept. In particular, in social sciences we 
have plenty of examples when people change their social 
roles, professions, status etci^. The scheme provides a 
straightforward way to computer simulations, where the 
rates identify the probabilities of change of a single agent. 

The scheme is applied to construct a model phase tran- 
sition, where one of two phases is chosen: the selfish 
strategy and the cooperative strategy. These options are 
discussed in therms of a crowd or a queue, where pedes- 
trians push the person before or stand still and move only 
when others move. We demonstrated, that at the points 
of the phase diagram where strategies are changed, the 
relaxation time of reaching the stationary phase goes to 
infinity. This effect, known in statistical mechanics as 
the critical slowing down, is basically a consequence of 
an increase of spatial correlation at the critical point. It 
is known to be present in numerical simulations, but it 
is reproduced also in mean field solutions. If a measure- 
ment of the related human behaviour is possible, this 
effect could be detected. 

As we know, neither mean field theory nor computer 
simulations is a decisive method to state if a phase tran- 
sition really exists. This is important in particular in 
the context of a one-dimensional Ising model with short 
range interaction. However, neither the thermodynamic 
limit nor the thermal equilibrium are met in social ex- 
periments, where systems are always finite and never sta- 
tionary. What kind of conclusions can be drawn, then. 



from the mean-field one-dimensional model? It seems 
reasonable to expect that if the states m± are close to 
be absorbing, they should appear in any finite system 
after a sufficiently long time. Then, the thermodynamic 
limit should leave them unchanged. The formal trick is 
therefore first to pass with time to infinity, later with the 
number of agents to infinity. However, even if the states 
m± are not fully absorbing, we should observe for appro- 
priate areas in the parameter space long-lasting periods 
when the majority of pedestrians plays this or that strat- 
egy. The situation is somewhat similar to the persistence 
of the stripes in the ground state of a two-dimensional 
Ising ferromagnet, which can be annealed out in a very 
small temperature^^. 

Concluding, all thermodynamic Monte Carlo simula- 
tions of magnetic phase transitions can be converted to 
athermal models, where the rates of switching of spins are 
treated as independent parameters. The rates are useful 
in Monte Carlo simulations, but they also enter to time- 
dependent mean field equations in a natural way. The 
athermal models can be used in areas other than physics 
of magnetism. For a sociologically motivated example of 
norm breaking in crowd, we have shown that the mean 
field approach reproduces the critical slowing down effect 
in the critical region. 
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